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Abstract. We develop further the approach to upper and lower bounds in 
quantum dynamics via complex analysis methods which was introduced by us 
in a sequence of earlier papers. Here we derive upper bounds for non-time av- 
eraged outside probabilities and moments of the position operator from lower 
bounds for transfer matrices at complex energies. Moreover, for the time- 
averaged transport exponents, we present improved lower bounds in the special 
case of the Fibonacci Hamiltonian. These bounds lead to an optimal descrip- 
tion of the time-averaged spreading rate of the fast part of the wavepacket in 
the large coupling limit. This provides the first example which demonstrates 
that the time-averaged spreading rates may exceed the upper box-counting 
dimension of the spectrum. 



1. Introduction 

This paper studies the long-time behavior of the solution to the time-dependent 
Schrodinger equation, idtip = Hip, in the Hilbert space -^^(Z) with an initially 
localized state, say ?/'(0) — Sq. More precisely, the Schrodinger operator is of the 
form 

(1) [Hu]{n) = u{n + 1) + u{n - 1) + V{n)u{n), 
where V : Z, ^ R is the potential, and the solution is given by 

(2) m = e-"«5o. 

The probability of finding the state in {n G Z : n > N} at time t is given by 

ri>Af 

Similarly, 

n<-N 

is equal to the probability of finding the state in {n G Z : n < —N} at time t. 

It is in general a hard problem to bound these so-called outside probabilities from 
above. Our recent paper |,13j introduced a way of estimating their time-averages 
from above in terms of the norms of transfer matrices at complex energies. In this 
paper we show how similar estimates can be obtained without the need to take 
time-averages. 
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To state these results, let us recall the definition of the transfer matrices. For 
71 G Z and z G C, define the transfer matrix <i>(n, z) by 

{T{n,z)- ■ ■T{l,z) n>l, 
Id n = 0, 

[T(n+l,z)]-i...[r(0,z)]-i n<-l, 

where 

z — V{m) —1 
1 

The definition is such that u : Z ^ C solves 



T(m, z) — 



(4) u{n + 1) + u{n - 1) + V{n)u{n) = zu{n) 

if and only if 



u{n) j ""^""'^M u(0) 



for every n £ Z. 



Theorem 1. Suppose H is given by ([T]), where V is a bounded real-valued function, 
and K > A is such that cr{H) C [—K + l,K — 1]. Then, the outside probabilities 
can be bounded from above in terms of transfer matrix norms as follows: 

i-K / \ -1 



Pr{N,t) <exp{-cN) + t* J ^^^max J|$(n,S + it~i)f^ dE, 
Pi(N,t) < cxp(-cN) + t^ [ ( max \\^ (n, E + it-^)\f^ dE. 



The implicit constants depend only on K . 

This result, which is the non-time averaged analogue of [T31 Theorem 7], has a 
number of consequences akin to those discussed in |13| . For p > 0, consider the 
p-th moment of the position operator, 

mi){t)^Y.\^n^^~''"^^^^r.)? 

and the upper and lower transport exponents (p) and (p) , given, respectively, 

/3,„(p) = hmsup^^^. 

and 

log(|X|^ ){t) 
/37(p) = liminf ^^-^ilil^. 

Let us briefly discuss a connection with the outside probabilities. Set 

P{N,t)^PiiN,t) + PriN,t). 
Following [1^, for < a < oo, define 

logP(r -l,t) 

t-^oo logt 

and 



(5) S (0;)=— liminf 

t — >oc 

and 

(6) 5'''' (a) = — lim sup 



logt 
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For every a, < 5+ (a) < S^{a) < oo. 

These numbers control the power decaying tails of the wavepacket. In particular, 
the following critical exponents are of interest: 

(7) = sup{a > : S'±(a) = 0}, 

(8) = sup{Q! > : S'=^(a) < oo}. 

We have that < < a~ < 1, < a^"^ < a+ < 1, and also that < , 
Q-u — • One can interpret af^ as the (lower and upper) rates of propagation 
of the essential part of the wavepacket, and a J as the rates of propagation of the 
fastest (polynomially small) part of the wavepacket; compare J16j. In particular, if 
a > a+, then P{t°',t) goes to faster than any inverse power of t. Since a ballistic 
upper bound holds in our case (for any potential V), a slight modification of [161 
Theorem 4.1] yields 



± 



lim Pf(p)^ai. 



and 

In particular, since (p) are nondecreasing, we have that 

(9) Pf^^ (p) < for every p>0. 

Corollary 1. Suppose H is given by ([ij, where V is a bounded real-valued function, 
and K > A is such that cr(H) C [—K + 1, if — 1] . Suppose that, for some C G (0, oo) 
and a G (0, 1), we have 

(10) / ( max \\'^(n,E + it-^)f'] dE = Oit-"") 



-K 



l<n<Ct° 



K / \ -1 

~1\ 



and 

(11) / ( max 11$ + ) dE ^ 0{t 
for every m> 1. Then 

(12) a+ < a. 
In particular, 

(13) f^to^P) — ^ f^^ every p > 0. 

Theorem [T] and Corollary [T] provide the analogues of the central general results 
from [13j for quantities that are not time-averaged. For comparison purposes, let 
us introduce the following notation. If f{t) is a function of i > and T > is 
given, we denote the time-averaged function at T by (/)(T): 

(/>(r) = |y^ e-''/^fit)dt. 

Thus, we consider, for example, {Pi{N, ■)){T), (P^(iV, •))(T), and {{\X\l)){T). For 
these time-averaged quantities, we can then define the transport exponents (/3J^^ (p)) 
and their limiting values (a^) and (a^) in the same way as above. For example, 
if in the formulation of Corollary [TJ we replace a+ by (a+) and /3^(p) by (/^^^(p)), 
we obtain the assertion of [13l Theorem 1]. 
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Let US now turn to a discussion of a special case. The Fibonacci Hamiltonian is 
the discrete one-dimensional Schrodinger operator in P (Z) as in ([1]) with potential 
: Z — > M is given by 

(14) V{n) = Ax[i-0-M)(n0~'+ mod 1). 

Here, A > is the coupling constant, (j) is the golden mean, 

and G [0, 1) is the phase. This is the most prominent model of a one-dimensional 
quasicrystal; compare the survey articles [UlJlj. It is known that the spectrum of 
H is independent of 6* let us denote it by Sa. Moreover, the Lebesgue measure 
of Ti\ is zero [23j and all spectral measures are purely singular continuous [S] . 

The quantum evolution with H given by the Fibonacci Hamiltonian has been 
studied in many papers. It had long been expected to be anomalous in the sense 
that it is markedly different from the behavior in the periodic case (leading to 
ballistic transport) and the random case (leading to dynamical localization); see, 
for example, papers in the physics literature by Abe and Hiramoto [T1[TB]. Lower 
bounds, showing in particular the absence of dynamical localization, were shown 
in [3l [71 [TOl [Til mi [m Ho] • There are far fewer paper establishing upper bounds 
for this model, especially for quantities like the moments of the position operator. 
Killip et al. showed for 9 — and A > 8 that the slow part of the wavepacket does 
not move ballistically [2^. Their result was extended to general in 5 . The first 
result establishing for A > 8 bounds from above for the whole wavepacket, and 
hence quantities like the moments of the position operator, is contained in jl3| . In 
that paper only the case 6* = is studied but it is remarked that the ideas from [5] 
will allow one to treat general 6''s. 

The paper ^13j introduced a new tool, the complex trace map, that allows one 
to use complex analysis methods in a context where real analysis methods were 
used earlier. This was important in our proof of non-trivial upper bounds for the 
(time-averaged) transport exponents. Given the analysis of 13 and Corollary [T] 
above, which strengthens [T3l Theorem 1], we obtain the following strengthening of 
[I3l Theorem 3]. 

Theorem 2. Consider the Fibonacci Hamiltonian, that is, the operator ([1]) with 
potential ()14|) . Assume that A > 8 and let 

2\og<t> 



a{\) 



\og Si{\y 



with 

(15) Si{\) = i ((A - 4) + v/(A-4)2-12) . 

Then, < a(A), and hence P^{p) < a(X) for every p > 0. 

As a byproduct of our study of the complex trace map in [13j , we established a 
distortion result that is useful to bound the transport exponents from either side. 
Since [13j focused on quantum dynamical upper bounds, we present the application 
of the distortion result to quantum dynamical lower bounds here. For this result, 
we still need to consider time-averaged quantities! 
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Theorem 3. Consider the Fibonacci Hamiltonian, that is, the operator ([T]) with 
potential (|14|). Suppose A > V24 and let 

SuiX) = 2A + 22. 

We have 

^''^ ^^^>^^^i^^is:ix)-pV'-i^^^ 

for a suitable constant C , and therefore 

21og0 



(17) «) > 



log5„(A)- 



Remark. In the special case = 0, the lower bound for {Pf^ip)) can be improved; 
see Theorem 0] at the end of this paper. 

The best previously known lower bound for (a,^) in the large coupling regime 
was obtained in [6J. It reads 

(18) «)>dim±(SA). 

Here, (Mu\^{T,\) denotes the upper/lower box counting dimension of Sa- Such 
a bound holds whenever the transfer matrices are polynomially bounded on the 
spectrum, which in particular holds in the Fibonacci case. 

The authors of [51 performed a detailed study of these dimensions. A particular 
consequence of their study is the following asymptotic statement, 

(19) lim dim±(EA) -logA^/^log^, 

A — ^oo 

where is an explicit constant (the unique maximum of an explicit function) that 
is roughly given by 

f* K 1.83156. 

On the other hand, for A > 8, [131 Theorem 3] (see also Theorem [2] above) gives 
the following upper bound for (aj), 

<"-> <- 

Thus, combining pT|) and (PO)) . we find the following exact asymptotic result. 
Corollary 2. For the Fibonacci Hamiltonian, we have 
(21) lim (aj) •logA = 21og0. 

A — >oo 

A particular consequence is that, as A — > oo, the limit behavior is the same for 
both (a^) and {a~). It would be of interest to show that these quantities are in 
fact equal for finite A. 

Comparing the asymptotic results (|19p and (HH), we see that the Fibonacci 
Hamiltonian at large coupling may serve as an example where the inequality in 
(HH]) is strict. This result is also relevant to a question raised by Last in [3T1 Sec- 
tion 9] who asked whether the upper box counting dimension of the spectrum could 
serve as an upper bound for dynamics and suggested that the expected answer is 
negative in general. See [171 125] for numerical results providing evidence supporting 
this expectation. 
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Let US put the recent results for the Fibonacci Hamihonian in perspective. The 
present paper and [6] result from an attempt to describe dimensions and transport 
exponents exactly. This is certainly a challenging problem and we have precise 
results only in an asymptotic regime; compare (jl9p and (j21[) . It would certainly 
be of interest to prove exact results also for fixed finite A. Moreover, the behavior 
in the small coupling regime is poorly understood. Many results that hold for A 
above some critical coupling do not obviously extend to smaller values since parts 
of the proofs break down. There is a definite need for new insights in order to prove 
results at small coupling. 

2. Upper Bounds for Outside Probabilities Without Time-Averaging 

Our goal in this section is to prove Theorem [T] and Corollary [TJ In particular, 
we will see how a specific consequence of the Dunford functional calculus allows us 
to replace the Parseval formula, which was a crucial ingredient in several earlier 
papers on quantum dynamical bounds [TOl [Til [HI [13 [HI [23 • This observation is 
the key to obtaining upper bounds for non-averaged quantities. 

Lemma 1. We have 

{e-''"So, <S„) = ^ e-^'^H ~ z)-^So, <5„) dz 

for every n S Z, i S M, and positively oriented simple closed contour T in C that is 
such that the spectrum of H lies inside F. 

Proof. This is a consequence of the so-called Dunford functional calculus; see [14] 
and also [H [22]. □ 



Lemma 2. Suppose H is given by ([T]), where V is a bounded real-valued function, 
and K>A is such that a{H) C [-~K + 1, if - 1]. Then, 

Pr{N,t)<exp{-cN)+ / ^ |((ff-S-ii-i)-Mo,(5„)pd£^, 

n>N 

Pi{N,t)<exp{-cN)+ J2 \{{H-E-it-Y'So,5,,)\''dE. 

Proof. Given i > 0, we will consider the following contour F: F = Fi U F2 U F3 U 
where 

Ti^{z = E + iy:Ee [-K,K], y = t-^} 

V2 = {z = E + iy:E=-K, y^[-l,t-^]) 

T:,^{z^E + iv:E^[-K,K], y = -1} 

Ti^{z^E + iy:E^K, ye [-1,^^]}. 

Notice that for 2; e F, we have ^z < t~^ and hence |e~'*^| = e*'^^ < e. Thus, by 
Lemma [TJ 

\{e'''"So,Sn)\<y2 I \{{H-z)-'5o,5n)\ \dz\. 

If z S F2 U F3 U F4, then clearly dist(z, cr(ii)) > 1 and hence the well-known 
Combes- Thomas estimate allows us to bound the contributions from F2 , F3 , F4 to 
PriN,t) and PiiN,t) by Cexp(-ciV). 
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The integral over Fi can be estimated using the Cauchy-Schwarz inequahty: 
(^^ \{{H~zr^do,S^)\\dz\^ <C{K) p^\{{H-E-it-^)-Ho,5^)\^dE. 

Combining these estimates, we obtain the assertion of the lemma. □ 
Proof of Theorem\^ We have the following estimates [13, p. 811]: 

(22) (P,(7V,.))(T)<exp(-ciV)+T-i /"^ V \{{H - E - iT-^)-Ho,5n)\^dE, 



>N 

(23) {Pi{N,-)){T)<eM-cN)+T-' r V \{{H - E - iT-^)-^do,6r^)\^dE. 
It was then shown in [13] how to derive 

{PriN,-)){T) <exp{-cN) J ^^^nmx J|$ (n, £; + iT^i) f ^ dE, 

{Pi{N,-)){T) <e^p{-cN) J ^ ^rnax ^J|$ (n, £; + iT-^) f ^ dE. 

from ([22|) and ((23|) . If we use Lemma [2] instead of ([22]) and ([23|) and then follow 
the very same steps, we obtain the desired bounds 

Pr{N,t)<eyip{-cN) + t^ j ^^^max J|$ (n, + rf-^) f ^ dE, 
PdN,t)<ex:n(-cN) + t^ ( ( max \\^ (n, E + it-^)\\A dE. 

J-K V-JV+1<»<0" ^ '"J 

This concludes the proof. □ 

Proof of CoroUaryUl Let us choose N{t) ~ [CPJ , where C E (0, oo) and a E (0, 1) 
are chosen such that ^ and ^ hold. Observe that P{N{t),t) = P{\N{t)~\,t). 
Then Theorem [1] shows that Pr{N{t),t) and Pi{N{t),t) go to faster than any 
inverse power of t. By definition of 5"+ (a) and (cf. ^ and (O), it follows that 
< a, which is (HI]). Finally, ([13]) follows from □ 

3. Tight Lower Bounds for Strongly Coupled Fibonacci 

3.1. The Complex Trace Map and the Distortion of Balls. In this subsection 
we recall some ideas from |13| and present an improvement of the key distortion 
result contained in that paper. 

For z e C and n g Z, consider the transfer matrices <I>(n, z) associated with the 
difference equation (JJ]) where V is the Fibonacci potential given by ([Ti]) . Notice 
that these matrices depend on both A and 6. 

Define the matrices Mk{z) by 

(24) ^g^oiFk,z) = Mk{z), fc>l, 

where Fk is the k-th Fibonacci number, that is, Fq = Fi = 1 and Fk+i = Fk + Fk-i 
for fc > 1. 

It is well-known that 

Mk+i{z)^ Mk-i{z)Mk{z), k>2. 



8 



D. DAMANIK AND S. TCHEREMCHANTSEV 



For the variables Xk{z) ~ iTrMfe(z)), fc > 1, we have the recursion 

(25) Xk+iiz) = 2xk{z)xk-iiz) - Xk-2(z) 
and the invariant 

(26) Xk+i{zf + Xkizf + Xk-i{zf - 2xk+i{z)xk{z)xk^i{z) - 1 = — . 

Letting x_i(z) = 1 and xq{z) = z/2, the recursion (I25p holds for all fc > 0. 
For ^ > 0, consider the sets 

fff = {z e C : \xk{z)\ < l + 

We have 

(^k U crf_i D crf+i U erf ^ Sa- 

Assume A > Ao(<5), where 

(27) Ao(^) - [12(1 + 5f + 8(1 + 5f + 4]i/2. 
The invariant (|26p implies 

(28) af n n af^^ = 0. 

Moreover, the set has exactly Fk connected components. Each of them is a 
topological disk that is symmetric about the real axis. 

It is known that all roots of a;^; are real. Consider such a root z, that is, Xk{z) = 
and define 

to(z) = #{0 < Z < fc - 1 : \xi{z)\ < 1}. 

Let 

Cfe,m = #{roots of Xk with m{z) — to}. 

An explicit formula for Ck.m was found in (BJ Lemma 5] (noting that our Ck,m 
equals ak^m + bk.m in the notation of that paper). All we need here is the following 
consequence of this result. 

Lemma 3. For every k > 2, Ck.m is non-zero if and only if ^ < m < 

Let {zl^^}i<j<Fk be the zeros of Xk and write to^-'^ — m{z]^'^) for I < j < Fk- 
Denote B{z,r) — {w G C : {w — z\ < r}. Then, we have the following distortion 
result. 

Proposition 1. Fix k > 3, 6 > 0, and A > max{Ao(2i5), 8}. Then, there are 
constants cs,ds > such that 

(29) |Ji?(z(^\ri'))cafc|jB(z(^),i?« 

where r'^^ = C55„(A)-"fc'\ and r'-/^ = ^^^((A)-™^'' . The first inclusion m ^ 
only needs the assumption A > Ao(2(5). 

Proof. The proof of this result is analogous to the proof of [TT, Proposition 3]. As 

explained there, TO^r*^ governs the size of |a;Jj(2:|,"''')|, which in turn is closely related to 

the size and shape of the connected component of erf that contains z^,"'-' by Koebe's 
Distortion Theorem. □ 
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3.2. Power-Law Upper Bounds for Transfer Matrices. In this subsection we 
prove power-law upper bounds for the norm of the transfer matrices for suitable 
complex energies and suitable maximal scales. 

Proposition 2. For every X,d > 0, there are constant C,"f such that for every k, 
every z £ cr^, and every N with 1 < N < Fk, we have 

(30) ||$(7V,z)|| < CAT^. 

The constant 7 behaves like O(logA) as A —> 00 for each fixed S. 

Proof. For 9 — 0, a modification of the proof of jljj Proposition 3.2] gives the result. 
The extension to arbitrary then follows using ideas from ^ . For the convenience 
of the reader, we sketch the main parts of the argument. 

Using the invariant ([26|) and [131 Lemma 4] , it is readily seen that 2 G implies 
the uniform bound 

(31) \xj{z)\<Cx,s forO<j<fc-l, 

where Ca,5 is an explicit constant that behaves like 0(A) as A 00; compare |1H 
Lemma 3.1]. 

Next, the uniform trace bound (|3ip yields an upper bound for the matrix norm, 

(32) l|M,(^)ll < {Cx.sY for < J < fc - 1, 

compare HT, Eq. (48)]. 

The final step is an interpolation of p2)) . Given any transfer matrix, ^{N, z), it 
is possible to write it as a product of matrices of type Mj{z); compare (SJ Section 3]. 
A careful analysis of the factors that occur, together with the estimate ([32]) , then 
gives (|30l) . The bound on 7 follows as in JTIJ Proposition 3.2]. □ 

3.3. Lower Bounds for the Transport Exponents. In this subsection we prove 
Theorem [3] using the results from the previous two sections. 

Proof of Theorem\^ By assumption, A > \/24. Thus, it is possible to choose 5 > Q 
so that A > Ao(2^). Recall that erf = {2; e C : |2;a;(-z)| < 1 + (5}. It follows from 
Lemma [3] and Proposition [1] that af, has a connected component Dk such that 

B{zk,rk) C Dk, 

where G M and 

rk = csSui\)-^^\ 
It follows from Proposition [2] that for z £ Dk and 1 < < F^, we have 

\\'^{N,z)\\ < CN-' 

with suitable constants C and 7. 

Let us fix k and take N = Fk (f)^ and e — ^ < Thus, 

r> — = -5„(A)r$i >iv^ 

rk cs 

where 

_ logg„(A) 
21og0 
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Due to the Parscval formula (see, e.g., [201 Lemma 3.2]), we can bound the 
time-averaged outside probabihties from below as follows, 

(33) {PiN,-)){T)>e [ \\^{N,E + te)\\-^dE. 

Jr 

To bound the integral from below, we integrate only over those E G (zk — rk, Zk+rk) 
for which E + ie G B{zk, r^) C D^. The length of such an interval 1^ is larger than 
crfc for some suitable c > 0. For E E 1^, we have 

\mN,E + ie)\\ <N'><Ti. 

Therefore, ([SS]) gives 



(34) (P(N,-)){T)>^T-^>T- 

where N = Fk, T > , for any k. 

Now let us take any T > 1 and choose k maximal with < T. Then, 

F^<T< < rF^. 

It follows from that 

{P{^Ti,-)){T) > {P{N, ■)){T) > 
for all T > 1. It follows from the definition of {j3^{p)) and {a~) that 



(/37ob)>>---0 



1 

and 

(a;> > -• 
s 

Since 

1 21og0 



s logs, (A)' 

and 7 = 0(log A), this completes the proof of Theorem[3l □ 

3.4. Improved Lower Bounds for Zero Phase. We conclude with a discussion 
of the special case 9 = Q. For this particular value of the phase, the result above can 
be improved. The key change is that in this case, one has "all the possible squares 
adjacent to the origin." This then enables one to directly estimate (P(iV, •))(T) 
from below without employing the general upper bounds for the transfer matrices. 

Theorem 4. Consider the operator ([1]) with potential (|14p and phase 6 ~ Q. Sup- 
pose A > ^/2A and let 

Su{\) = 2A + 22. 

We have 

2\og4> 2 /', 21og4I \ 



for every p > 0. 
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Proof. Starting again with the Parscval formula [201 Lemma 3.2], we have 

(36) (P,(iV,.))(T)^-l ^ / \{{H-E~iT-Y^6o,S.r.)\'dE. 

n>N''^ 

Using the fact that 

u(n) - {{H -E-iT-^5Q,5n) 

solves the difference equation Hu = {E+iT~^)u away from zero, one can bound the 
right-hand side of (j36p from below by the Gordon-type mass reproduction technique 
used, for example, in [3] [71 [121 [19] . We refer the reader to [l2l Lemmas 2 and 3] 
for a formulation of this technique that is particularly well suited to our situation 
here. It is important to point out that while in those papers, only real energies are 
considered, the statements extend to complex energies as soon as one can arrange 
for the same transfer matrix trace estimates to hold. 

As in the proof of Theorem [H since A > \/24, it is possible to choose (5 > so 
that A > Xo{2S). Fixing this value of i5, we consider the set af, — {z £ C : \xk{z)\ < 
1 4-i5}. Notice that Xk{z) is one-half the trace of the transfer matrix of the potential 
under consideration from the origin to Fk. This is why the case = is special! It 
follows from [13l Lemma 4] that for < j < fc — 1, min{|a;j(2:)|, |a;j+i(z)|} < 1 + S 
for every z af.. Since, for every j, the potential repeats its values from 1 to Fj 
once, we can apply [12l Lemma 2] to bound the right-hand side of from below. 
Of course we use either j or j -I- 1, depending on which of \xj{z) and |a;j_|-i(z)| is 
bounded by 1 + (5. Carrying this out inductively, we find that 

^ \uin)f>CN'^i\u{l)\' + \uiO)f) 

n>N 

for a (5-dependent constant 

K < — 

51og0 

and a corresponding suitable constant C > 0. One checks that k can be made 

log 

arbitrarily close to ^j^^"^ by making S sufficiently small. The same reasoning can 
be applied to {Pi{N, ■)){T). In the end, one only needs to observe that 

\uii)\^ + \um' + i«(-i)p > c(i + i^F{E + ±)n 

uniformly in the energy, where F is the Borel transform of the spectral measure; 
compare [Till [TT] . 

From this point on, one can proceed as in the proof of Theorem [3l It then 
follows, with the same notation as above, that 

(37) (P(cT^,-))(T) >cT-i+2'=/^ / {l + {QF(E+^)f)dE 
and hence 

mp))>---(^--) 

° s p \ sJ 

from which the result follows. □ 
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Remark. It is possible to improve the estimate psp if one revisits the inequahty 
P7p and observes that the spectral measure is singular and the integral over Ik in 
question is larger than just This can be made quantitative using ideas from 
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